AD-A148  784  NUMERICAL  COMPUTATION  OF  BASE  FLOW  FOR  A  MISSILE  IN  THE 
PRESENCE  OF  A  CEN.  .  (U)  ARMV  BALLISTIC  RESEARCH  LAB 
ABERDEEN  PROVING  GROUND  MD  J  SAHU  ET  AL.  OCT  84 
UNCLASSIFIED  BRL-NR-3297  SBI-AD-F288  518  F/G  2074 


- 

JL 

A 

r 

- 

* 

- 

r. 

MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS -1963 -A 


AD- A 148  784 


AD£  s/o 


B 

R 

L 


MEMORANDUM  REPORT  BRL-MR-3397 


j 


NUMERICAL  COMPUTATION  OF  BASE  FLOW 
FOR  A  MISSILE  IN  THE  PRESENCE 
OF  A  CENTERED  JET 


f'M,  Eu 


‘J  * 
*  . 


Jubaraj  Sahu 
Charles  J.  Nietubicz 


October  1984 


NOV  2  7 1984 


APPROVED  FOR  PUMJC  RELEASE;  DISTRIBUTION  UNLIMITED. 

\  era 


US  ARMY  BALLISTIC  RESEARCH  LABORATORY 

ABERDEEN  PROVING  GROUND,  MARYLAND 


84  11  26  243 


Destroy  this  report  when  it  is  no  longer  needed. 
Do  not  return  it  to  the  originator. 


Additional  copies  of  this  report  may  be  obtained 
from  the  National  Technical  Information  Service, 

U.  S.  Department  of  Commerce,  Springfield,  Virginia 


22161. 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position,  unless  so  designated  by  other 
authorized  documents. 


The  use  of  trade  names  or  manufacturers'  names  in  this  report 
does  not  constitute  indorsement  of  any  commercial  product. 


SECURITY  CLASSIFICATION  OF  THIS  RACE  r«hai  Dim  entered) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


rmwiriTin 


MEMORANDUM  REPORT  BRL-MR-3397 


4.  TITLE  rand  Submit 


NUMERICAL  COMPUTATION  OF  BASE  FLOW  FOR  A  MISSILE 
IN  THE  PRESENCE  OF  A  CENTERED  JET 


7.  AUTHONfA) 

Jubaraj  Sahu  and  Charles  J.  Nietubicz 


S.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

U.S.  Arny  Ballistic  Research  Laboratory 
ATTN:  AMXBR-LFD 

Aberdeen  Proving  Ground,  Maryland  21005-5066 


It.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

U.S.  Army  Ballistic  Research  Laboratory 
ATTN:  AMXBR-OD-ST  .  , 

Aberdeen  Proving  Ground.  MD  21005-5066 


S.  TYPE  OF  REPORT  A  PERIOO  COVERED 


Final 


S.  PERFORMING  ORG.  REPORT  NUMBER 


S.  CONTRACT  OR  GRANT  NUMBER!*; 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  4  WORK  UNIT  NUMBERS 


RDT&E  1L161102AH43 


12.  REPORT  DATE 

October  1S84 


is.  number  of  pages 


MONITORING  AGENCY  NAME  *  ADORESSTI/  dlttorant  hem  Controllln*  Olllco)  I  IS.  SECURITY  CLASS,  (ol  O Urn  report) 


Unclassified 


IS.  DISTRIBUTION  STATEMENT  (ol  tflta  R* port) 


Approved  for  public  release,  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ol  the  obattoct  entered  In  Block  30,  II  different  hem  Report) 


It.  SUPPLEMENTARY  NOTES 


This  report  supersedes  BRL  IMR  789,  dated  September  1983. 


IS.  KEY  WOROS  (Continue  on  TO  retro  ride  II  nocoooory  and  Identify  by  Slock  number; 

Navler- Stokes  Computations  Base  Flow 

Propulsive  Jet  Compressible 

Implicit  Algorithm  Time  Dependent 

Axi symmetric  Afterbody 


I  »0>  A— TWACT  CCmtlmmmm  mw 


*  mmd  Mmtiitr  Wf  block  numbm) 


A  thin-layer  Navier-Stokes  base  flow  code,  developed  for  projectile  aero¬ 
dynamics,  has  been  used  to  compute  the  missile  base  flow  in  the  presence  of  a 
centered  propulsive  jet.  The  thin-layer,  compressible,  Navier-Stokes 
equations  are  solved  using  a  time  dependent,  implicit  numerical  algorithm. 
Numerical  computations  have  been  made  to  predict  the  missile  afterbody  flow 
field  for  both  jet-off  and  jet-on  conditions.  Solutions  are  obtained  for  an 
axi symmetric  cylindrical  afterbody  where  the  free  stream  Mach  number  is  1.343 


FORM 

I  JAM  7* 


EDITION  OF  •  MOV  BB  IB  I 


UNCLASSIFIED 

iecurity  classification  of  this  page  fB 


i  Data  Intaratf) 


UNCLASSIFIED 


MCUMTV  CLMWFjCATjOW  OF  THU  PAOVUum  Of 

20.  ABSTRACT  (Continued) 

and  the  jet  exit  Mach  number  is  2.7.  The  exhaust  jet  static  pressure  is  2.15 
times  the  free  stream  static  pressure.  Computed  results  show  both  qualitative 
and  quantitative  features  of  the  base  region  flow  field. 


UNCLASSIFIED _ 

sccuftiTY  classification  of  this  PAGeriw»#n  Bnffd) 


v.  \ . 


•*.»*,  «*.  i*.  *’•  •*.  * 


TABLE  OF  CONTENTS 


I. 

II. 

II. 


IV. 

V. 

VI. 

II. 


LIST  OF  ILLUSTRATIONS . 

INTRODUCTION . 

COMPUTATIONAL  TECHNIQUE . 

METHOO  OF  SOLUTION . 

A.  Base  Region  Flow  with  Jet-Off 

B.  Base  Region  Flow  with  Jet-On. 

MODEL  AND  TEST  CONDITIONS . 

COMPUTATIONAL  GRID . 

RESULTS . 

SUMMARY . 

REFERENCES . 

LIST  OF  SYMBOLS . 


DISTRIBUTION  LIST 


LIST  OF  ILLUSTRATIONS 

Figure  Page 

1  Schematic  Illustration  of  Flow  Field  Segmentation .  16 

2  MICOM  Base-Flow  Model  Mounted  In  the  Wind  Tunnel .  17 

3  Base  Flow  Field  Measurement  Regions .  18 

4  Configuration  and  Notation  for  the  Experimentally  Determined 

Boundary  Data . . .  19 

5  Computational  Grid .  20 

6  Expanded  View  of  the  Grid  Near  the  Model .  21 

7  Velocity  Vectors,  M„  «  1.343,  a  «  0,  Jet-Off .  22 

8  Velocity  Vectors,  Mw  *  1.343,  a  »  0,  Jet-On .  23 

9  Velocity  Vectors  with  Model  Drawn  to  Scale,  M„  =  1.343,  a  =  0, 

Jet-On .  24 

10  Stream  Function  Contours,  M„  *  1.343,  o  =  0,  Jet-Off .  25 

11  Stream  Function  Contours,  M^  *  1.343,  a  «  0,  Jet-On .  26 

12  Variation  of  Base  Pressure  Along  the  Base,  M  *  1.343,  a  =  0, 

Jet-Off .  27 

13a  Cone  Wind  Tunnel  Model . 28 

13b  Base  Pressure  Distribution  for  Cone .  29 

14  Expanded  Grid  for  the  Isolated  Base  Region  Flow  Field 

Computations . 30 

15  Variation  of  Base  Pressure  along  the  Base,  M^  =  1.343,  a  =  0, 

Jet-Off  and  Jet-On .  31 

16  Static  Pressure  Distribution  at  Z/D  =  .5,  M,,  =  1.343,  a  =  0, 

Jet-On . 32 

17  Pressure  Contours,  Mw  *  1.343,  a  *  0,  Jet-On .  33 

18  Mach  Contours,  Mm  *  1.343,  a  -  0,  Jet-On .  34 


I.  INTRODUCTION 


The  aerodynamic  drag  associated  with  afterbody  configurations  is  influ¬ 
enced  both  by  the  afterbody  geometry  and  by  the  interaction  of  the  exhaust  jet 
with  the  external  flow.  The  exhaust  jet  has  a  large  influence  on  the  base 
pressure  and  can  influence  the  flow  on  the  afterbody  for  large  jet  pressure. 
The  flow  field  for  such  cases  is  complex  due  to  strong  viscous/inviscid  inter¬ 
action  and  regions  of  flow  recirculation.  Because  of  the  strong  interaction 
involved  between  flow  regimes.  It  Is  advantageous  to  use  the  Navier- Stokes 
computational  technique  since  It  considers  the  Interactions  in  a  fully-coupled 
manner. 

As  part  of  a  continued  research  effort  at  BRL,  numerical  computational 
capabilities  have  been  developed  to  predict  the  aerodynamic  behavior  of  artil¬ 
lery  shell.  References  1  and  2  have  reported  the  development  and  application 
of  the  Azimuthal-Invariant  Thin-Layer  Navier-Stokes  computational  technique  to 
predict  the  flow  about  slender  bodies  of  revolution  at  transonic  speeds. 
These  calculations,  however,  did  not  Include  the  base  region  flow  field  and 
modeled  the  wake  as  an  extension  of  the  afterbody.  Recently,  a  new  numerical 
capability3’4  has  been  developed  for  the  simulation  of  the  complete  projectile 
including  the  base  region.  This  code  has  been  used  to  predict  the  base  pres¬ 
sure  of  shell  at  transonic  speeds  both  with  and  without  base  bleed.  The 
objective  of  the  present  research  work  was  to  use  this  base  flow  code  and 
extend  Its  capability  to  predict  the  effects  of  a  centered  propulsive  jet  on 
the  base  region  flow  field.  The  same  numerical  technique  has  been  used  by 
Deiwert5  for  a  computational  Investigation  of  supersonic  axisymmetric  after¬ 
bodies  containing  a  centered  propulsive  jet.  However,  his  method  of  solution 


1 .  C.  J.  Nietubicz ,  T.  H.  Pulliam,  and  J .  L.  Stager,  "Numerical  Solution  of 
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Projectile  Base  Flow  with  and  without  Base  Injection ,"  U.S.  Army  Ballistic 
Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland,  ARBRL-TR-02532, 
November  1983.  (AD  A133684) 

5.  G.  S.  Deiwert,  "A  Computational  Investigation  of  Supersonic  Axisymmetric 
Flow  Over  Boattails  Containing  a  Centered  Propulsive  Jet,"  AIAA  Paper  No. 
83-0462,  January  1983. 


and  treatment  of  the  computational  mesh  differs  substantially  from  our  current 
approach. 

A  brief  description  of  the  computational  technique  is  given  in  Section 
II.  A  novel  flow  field  segmentation  procedure,  equivalent  to  using  multiple 
adjoining  grids.  Is  used  and  described  in  Section  III.  The  model  geometry  and 
test  conditions  are  given  In  Section  IV  and  the  generation  of  the  computation¬ 
al  grid  is  described  in  Section  V.  In  Section  VI,  results  are  presented  for 
afterbody  flow  field  computations  for  both  jet-off  and  jet-on  conditions.  The 
free  stream  Mach  number  is  1.343  and  the  jet  exit  Mach  number  is  2.7.  The 
exhaust  jet  static  pressure  is  2.15  times  thefree  stream  pressure  and  the  con¬ 
ical  nozzle  half  angle  is  10°.  The  Reynolds  number  of  the  external  flow  based 
on  maximum  body  diameter  is  0.92  x  10®.  These  computations  correspond  to  the 
experiments  conducted  by  Walker6  for  the  same  configuration  and  flow  condi¬ 
tions.  The  fluid  has  been  assumed  to  be  air  in  both  the  external  flow  and  the 
exhaust,  although  nitrogen  was  used  in  the  exhaust  in  the  actual  experiments. 


II.  COMPUTATIONAL  TECHNIQUE 

The  Azimuthal  Invariant  (or  Generalized  Axi symmetric)  thin-layer  Navier- 
Stokes  equations  for  general  spatial  coordinates  n,  t  can  be  written  as2 

aTq  +  3eE  +  3;G  +  H  =  Re'^S  (1) 


where 


5  ■  £(x,y,z,t)  is  the  longitudinal  coordinate 
n  =  n(y,z,t)  is  the  circumferential  coordinate 
C  =  c(x,y,z,t)  is  the  near  normal  coordinate 
t  *  t  is  the  time 
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6.  B.  J.  Walker,  "Tactical  Mieeion  Base  Flaw,  "  Paper  presented  at  Symposium 
<xi  Rocket/Plume  Fluid  Dynamic  Interactions,  Huntsville,  Alabama,  April 
1983. 
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The  velocities 

U  =  5t  +  5XU  +  Syv  +  Czw 

V  =  nf  +  nu  +  n  v  +  n^w  (2) 

v  X  Jr  a 

W  =  ct  +  cxu  +  Cyv  +  czw 

represent  the  contravariant  velocity  components. 

The  Cartesian  velocity  components  (u,  v,  w)  are  nondimensional ized  with 
respect  to  a,,  (free  stream  speed  of  sound).  The  density  (p)  is  referenced 
to  pm  and  total  energy  (e)  to  p^a2.  The  local  pressure  is  determined  using 
the  equation  of  state, 

P  *  (y  -  1 ) [e  -  0.5p(u2  +  v2  +  w2)]  (3) 

where  y  is  the  ratio  of  specific  heats. 

In  Equation  (1),  axisymmetrlc  flow  assumptions  have  been  made  which 

A 

result  in  the  source  term,  H.  The  details  of  how  this  is  obtained  can  be 
found  in  Reference  2  and  are  not  discussed  here.  Equation  (1)  contains  only 
two  spatial  derivatives.  However,  it  retains  all  three  momentum  equations  and 


allows  a  degree  of  generality  over  the  standard  axl symmetric  equations.  In 
particular,  the  circumferential  velocity  is  not  assumed  to  be  zero^thus  allow¬ 
ing  computations  for  spinning  projectiles  to  be  accomplished. 

The  numerical  algorithm  used  Is  the  Beam-Warming  fully  implicit,  approxi¬ 
mately  factored  finite  difference  scheme.  The  algorithm  can  be  first  or 
second  order  accurate  In  time  and  second  or  fourth  order  accurate  in  space. 
Since  the  Interest  Is  only  In  the  steady-state  solution.  Equation  (1)  is 
solved  In  a  time  asympotic  fashion  and  first  order  accurate  time  differencing 
Is  used.  The  spatial  accuracy  Is  fourth  order.  Details  of  the  algorithm  are 
Included  In  References  7-9. 

For  the  computation  of  turbulent  flows,  a  turbulence  model  must  be 
supplied.  In  the  present  calculations  a  two  layer  algebraic  eddy  viscosity 
model  by  Baldwin  and  Lomax10  is  used.  In  their  two  layer  model  the  inner 
region  follows  the  Prandtl-Van  Driest  formulation.  Their  outer  formulation 
can  be  used  In  wakes  as  well  as  In  attached  and  separated  boundary  layers.  In 
both  the  inner  and  outer  formulations,  the  distribution  of  vorticity  is  used  to 
determine  length  scales,  thereby  avoiding  the  necessity  of  finding  the  outer 
edge  of  the  boundary  layer  (or  wake).  The  magnitude  of  the  local  vorticity 
for  the  axl symmetric  formulation  Is  given  by 


M  *  /  <£>2  +  « -  W  *  « -  £>2 . 


(4) 


In  determining  the  outer  length  scale  a  function10 


F(y)  -  y|“|  [1  -  exp(-y+/A+)] 


(5) 


7.  J.  L.  Steger,  "Implicit  Finite  Difference  Simulation  of  Flow  About 
Arbitrary  Geometries  with  Application  to  Airfoils AIAA  Journal .  Vol. 

16,  No.  7,  July  1978,  pp.  679-686. 

8.  T.  H.  Pulliam,  and  J.  L.  Steger,  "On  Implicit  Finite-Difference  Simula¬ 
tions  of  Three-Dimensional  Flaw, "  AIAA  Journal.  Vol.  18,  No.  2,  February 
1980,  pp.  159-167. 

9.  R.  Beam,  and  R.  F.  Warming,  "An  Implicit  Factored  Scheme  for  the  Compres¬ 
sible  Navier-Stokes  Equations."  AIAA  Journal.  Vol.  16,  No.  4,  April  1978, 
pp .  393— 402 . 


10.  B.  S.  Baldwin,  and  H.  Lomax,  "Thin-Layer  Approximation  and  Algebraic 
Model  for  Separated  Turbulent  Flows,"  AIAA  Paper  No.  78-257,  1978. 


Is  used  where  y+  and  A+  are  the  conventional  boundary  layer  terms.  For  the 
base  flow  (or  wake  flow)  the  exponential  term  of  Equation  (5)  is  set  equal  to 
zero. 


III.  METHOD  OF  SOLUTION 
A.  Base  Region  Flow  with  Jet-Off 

The  procedure  used  to  compute  t..j  base  flow  without  jet  for  projectile 
configuration  has  been  described  in  Reference  3;  however,  limited  details  will 
be  repeated  here  for  clarity.  The  code  computes  the  full  flow  field  (includ¬ 
ing  the  base  region)  of  a  projectile.  Figure  1  shows  a  schematic  illustration 
of  the  flow  field  segmentation  used  in  this  study  for  computational  purposes. 
It  shows  the  transformation  of  the  physical  domain  into  the  computational 
domain  and  the  details  of  the  flow  field  segmentation  procedure  in  both 
domains.  This  flow  field  segmentation  procedure  is  equivalent  to  using  multi¬ 
ple  adjoining  grids.  An  important  advantage  of  this  procedure  lies  in  the 

preservation  of  the  sharp  corner  at  the  base  and  allows  easy  blending  of  the 

computational  meshes  between  the  regions  ABCD  and  AEFG.  No  approximation  of 
the  actual  sharp  corner  at  the  base  is  made.  Thus,  realistic  representation 
of  the  base  is  inherent  in  the  current  procedure. 

The  cross  hatched  region  represents  the  model.  The  line  BC  is  the  base 
and  the  region  ABCD  is  the  base  region  or  the  wake.  The  line  AB  is  a  computa¬ 
tional  cut  through  the  physical  wake  region  which  acts  as  a  repetitive  bound¬ 
ary  in  the  computational  domain.  Implicit  integration  is  carried  out  in  both 
£  and  c  directions  (see  Figure  1).  Note  the  presence  of  the  lines  BC  (base) 
and  EF  (nose  axis)  in  the  computational  domain.  They  both  act  as  boundaries 
in  the  computational  domain  and  special  care  must  be  taken  in  inverting  the 
block  tridiagonal  matrix  in  the  5  direction.  The  details  of  these  can  be 

found  in  References  3  and  4  and  are  not  included  here. 

The  no  slip  boundary  condition  for  viscous  flow  is  enforced  by  setting 


U  =  V  =  W  =  0  (6) 


on  the  body  surface.  At  the  base  boundary,  inviscid  boundary  condition  has 
been  used.  Additionally,  at  the  corner  of  the  base,  the  boundary  conditions 
are  double-valued  and  depend  on  the  direction  from  which  the  corner  is 
approached.  Approaching  it  in  the  streamwise  direction  (EB),  the  no-slip 
boundary  condition  is  used;  while  approaching  in  the  radial  direction  (along 
the  base,  CB),  inviscid  boundary  condition  is  used.  Along  the  computational 
cut  (AB) ,  the  flow  variables  above  and  below  the  cut  are  simply  averaged  to 
determine  the  boundary  conditions  on  the  cut.  On  the  centerline  of  the  wake 
region,  a  symmetry  condition  is  imposed  and  free  stream  conditions  are  used  on 
the  outer  boundary.  First  order  extrapolation  for  all  flow  variables  is  used 
at  the  downstream  boundary. 
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B.  Base  Region  Flow  with  Jet-On 


The  method  of  solution  for  the  case  with  a  centered  propulsive  jet  remains 
essentially  the  same  as  described  in  Section  III  A.  The  boundary  conditions 
on  the  body  surface,  the  cut, and  the  downstream  boundary  also  remain  the  same 
as  previously  described.  Along  that  part  of  the  base  boundary  which  does  not 
contain  the  nozzle,  same  conditions  described  earlier  are  used.  For  the 
remaining  part  of  the  base  boundary  (i.e. ,  the  nozzle  exit),  boundary  condi¬ 
tions  are  used  based  on  the  nozzle  exit  Mach  number,  stagnation  temperature, 
and  pressure.  The  velocity  components  are  linearly  interpolated  from  the 
center  line  of  symmetry  to  the  nozzle  height  at  the  exit  i.e.,  conical  flow  at 
the  jet  exit  has  been  assumed. 


IV.  MODEL  AND  TEST  CONDITIONS 

The  model  geometry  used  in  this  study  is  shown  in  Figure  2.  It  consists 
of  a  4  caliber  tangent-ogive  nose  and  a  5  caliber  cylindrical  afterbody. 

The  base  diameter  is  1  caliber  and  the  nozzle  exit  diameter  is  0.2  caliber. 

The  experimental  model  was  side-wall  mounted  in  AEDC  Wind  Tunnel6  1-T  as 
shown  in  Figure  2a;  additional  details  of  the  nozzle  and  base  are  shown  in 
Figure  2b.  The  wind  tunnel  and  nozzle  test  conditions  are  summarized  below. 


Free  Stream  Flow 


Gas:  Air 

Stagnation  Temperature:  190°F 
Static  Pressure:  6.14  psia 
Free  Stream  Mach  Number:  1.343 


Nozzle  Flow 
Gas:  Nitrogen 

Stagnation  Temperature:  207 °F 
Stagnation  Pressure:  307.7  psia 
Exit  Design  Mach  Number:  2.7 


Detailed  experimental  measurements  of  the  base  region  flow  field  have  been 
made  by  Walker.6  The  data  were  obtained  using  a  two-component  LDV  system. 
Figure  3  is  a  sketch  of  the  aft  end  of  the  model  and  shows  the  regions  for 
which  the  data  is  presented.  The  test  results  have  not  been  released.  A 
"sealed  envelope"  comparison  between  computational  models  and  experimental 
results  for  the  MICOM  model  with  a  centered  propulsive  jet  is  being  carried 
out  by  the  U.S.  Army  Missile  Command.  Various  groups  are  performing  computa¬ 
tions  and  MICOM  will  make  the  comparison.  At  this  point,  limited  experimental 
data  is  available  in  the  form  of  velocity  components  and  Mach  number  at  two 
boundary  stations  (see  Figure  4).  The  available  experimental  data  at  a 
station  1.6  caliber  upstream  of  the  base  has  been  used  as  an  upstream  boundary 
condition  In  some  of  the  numerical  calculations.  The  fluid  has  been  assumed 
to  be  air  for  bothfree  stream  and  nozzle  flows  in  the  present  computations. 


V.  COMPUTATIONAL  GRID 


Figure  5  shows  the  computational  grid  used  for  the  numerical  computations. 
This  grid  was  obtained  In  two  segments  consistent  with  the  flow  field  segmen¬ 
tation  procedure  described  earlier  in  Section  III  A.  The  first  segment  of  the 
computational  mesh  Is  for  region  AEFG  (see  Figure  1)  and  was  obtained  from  a 
grid  generation  program  developed  by  Steger  et  al.11  The  grid  for  the  second 
segment,  region  ABCD,  was  obtained  by  another  procedure  to  be  described  later 
In  this  section.  The  grid  generation  program  by  Steger  et  a?11  allows  arbi¬ 
trary  grid  point  clustering,  thus  enabling  grid  points  for  the  projectiles  to 
be  clustered  in  the  vicinity  of  the  body  surface.  The  grid  consists  of  109 
points  in  the  longitudinal  direction  with  25  points  in  the  base  region  and  40 
points  In  the  radial  direction.  The  computational  domain  extended  to  little 
more  than  a  body  length  in  front,  about  3  body  lengths  in  the  radial  direction 
and  about  2  body  lengths  behind  the  base  of  the  missile.  The  grid  points  in 
the  normal  direction  were  exponentially  stretched  away  from  the  surface  with  a 
minimum  spacing  at  the  wall  of  .00002D.  This  spacing  locates  at  least  two 
points  within  the  laminar  sublayer. 

An  expanded  view  of  the  grid  in  the  vicinity  of  the  model  is  shown  in 
Figure  6.  As  stated  earlier,  the  grid  shown  in  Figures  5  and  6  was  generated 
in  two  segments.  First,  the  grid  in  the  outer  region  AEFG  is  obtained  using 
an  elliptic  solver11  for  the  ogive  portion.  Straight-line  rays  are  used  for 
the  remaining  portion  which  runs  all  the  way  to  the  downstream  boundary. 
Second,  the  grid  in  the  base  region  ABCD  is  obtained  simply  by  extending  the 
straight  lines  perpendicular  to  line  AB  down  to  the  line  of  symmetry  (line 
CD).  In  other  words,  same  grid  stretching  is  used  in  the  longitudinal  direc¬ 
tion  in  both  the  regions  above  and  below  the  line  AB.  It  remains  then  to 
generate  the  grid  spaclngs  in  the  normal  direction  in  the  base  region  ABCD. 
For  base  flow  computations  with  a  centered  propulsive  jet,  it  is  desirable  to 
cluster  grid  points  near  the  cut  to  resolve  the  flow  gradients  in  the  shear 
layer  as  well  as  near  the  line  of  symmetry  where  jet  is  present.  A  1-D  ellip¬ 
tic  equation  solver  was  used  to  generate  the  grid  spacings  in  the  normal 
direction  in  the  base  region  ABCD.  This  technique  enables  the  desired  minimum 
spacing  to  be  specified  at  both  ends.  It  should  be  noted  that  the  same  mini¬ 
mum  spacing  .00002D  is  specified  on  both  sides  of  the  cut  thus  maintaining  a 
smooth  variation  of  grid  across  the  cut.  This  spacing  could,  of  course,  be 
increased  downstream  of  the  base.  The  number  of  grid  points  above  and  below 
line  AB  is  the  same  (40  points).  As  can  be  seen  in  Figure  6,  the  grid  points 
are  clustered  near  the  nose-cylinder  junction  and  at  the  projectile  base  where 
appreciable  changes  in  flow  variables  are  expected. 

When  computations  over  the  entire  model  are  made,  only  a  limited  number  of 
grid  points  can  be  used  in  the  base  region.  One  way  to  eliminate  this 
restriction  is  to  use  known  data  given  by  experiment  or  otherwise  at  a  station 
upstream  of  the  base  and  then  compute  the  flow  field  in  the  isolated  base 


11.  J.  L.  Steger,  C.  J.  Nietubicz,  and  K.  R.  Reavey,  "A  General  Curvilinear 
Grid  Generaticn  Program  for  Projectile  Configurations, "  U.S.  Army 
Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland,  ARBRL- 
MR-03142,  October  1981.  (AD  A107334) 
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region  only.  This,  of  course,  allows  a  large  number  of  grid  points  to  be  used 
In  the  base  region  and  can  be  used  to  determine  grid  dependancy  on  the  compu¬ 
ted  solution  in  the  base  region.  This  capability  was  initially  developed  by 
the  current  authors  for  the  numerical  computations  of  base  region  flow  field 
at  supersonic  velocities.  Solutions  were  obtained  for  a  cone-cylinder  config¬ 
uration  with  the  space-marching  PNS  (Parabolized  Navier- Stokes)  code.  This 
generated  a  solution  at  a  station  upstream  of  the  base,  and  this  solution  was 
then  used  as  the  upstream  boundary  condition  for  the  computation  of  the  base 
region  flow  field  by  the  unsteady  base  flow  code.  In  the  present  computa¬ 
tions,  experimental  data  is  known  at  a  station  upstream  of  the  base,  and  is 
thus  used  as  the  upstream  boundary  condition  in  the  numerical  computations 
when  full  grid  is  not  used. 


VI.  RESULTS 

Results  are  now  presented  for  both  jet-off  and  jet-on  cases.  Figure  7 
shows  the  velocity  vectors  in  the  base  region  when  there  is  no  jet  present. 
Each  vector  shows  the  magnitude  and  direction  of  the  velocity  at  that  point. 
The  recirculatory  flow  in  the  base  region  is  clearly  evident.  Figure  8  shows 
the  effect  of  the  centered  jet  on  the  near  wake  flow  field.  The  near  wake 
flow  field  has  changed  considerably  and  the  recirculation  region  has  been 
reduced.  This  figure  has  been  stretched  in  the  vertical  direction  to  show  the 
qualitative  features  as  clearly  as  possible.  A  scaled  result  is  shown  in 
Figure  9,  and  Is  a  more  realistic  representation  of  the  near  wake  flow  field 
with  jet-on. 

Figures  10  and  11  show  the  stream  function  contours  in  the  base  region. 
Figure  10  Is  for  the  jet-off  case  and  clearly  shows  the  separation  bubble. 
The  reattachment  point  is  about  two  calibers  down  from  the  base.  The  effect 
of  the  jet  is  shown  in  Figure  11.  The  separation  bubble  has  been  considerably 
reduced  in  size. 

Quantitatively,  one  is  Interested  in  how  the  complex  flow  field  in  the 
base  region  affects  the  base  pressure.  Figure  12  shows  the  variation  of  base 
pressure  along  the  base  for  jet-off  conditions.  When  there  is  no  jet  present, 
base  pressure  is  high  near  the  corner  of  the  base.  It  decreases  gradually  and 
then  recovers  to  almost  the  free  stream  static  pressure  near  the  center  line 
of  symmetry.  This  variation  of  base  pressure  along  the  base  is  substantiated 
by  other  experimental  observations.1*  The  cone  model  for  the  experimental 
measurements  of  Reference  12  is  shown  in  Figure  13a.  The  model  was  supported 
by  a  sting  attached  to  the  base.  Figure  13b  shows  the  typical  base  pressure 
variation  along  the  base  for  various  Reynolds  numbers  and  =  1.75.  There  is 
some  effect  of  the  sting  on  these  variations;  however,  the  experimental  base 
pressure  profiles  have  the  same  qualitative  variations  as  predicted  by  the 
present  computations  (Figure  12). 


12.  L.  D.  Kayser,  "Experimental  Study  of  Separation  from  the  Baee  of  a  Cone 
at  Supersonic  Speeds,"  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen 
Proving  Ground,  Maryland,  BRL  Report  No.  1737,  August  1974.  (AD  A005015) 
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The  primary  objective  of  this  research  was  to  compute  the  base  flow  in  the 
presence  of  a  centered  propulsive  jet.  For  the  jet-on  computations,  the  full 
grid  shown  In  Figure  6  was  used.  Additionally,  to  determine  grid  dependancy 
on  the  solution,  a  large  number  of  grid  points  were  used  in  the  base  region 
and  solutions  were  obtained  for  another  grid  shown  in  Figure  14.  In  this 
case,  calculations  were  not  performed  on  the  full  model.  The  given  experi¬ 
mental  data  was  used  as  an  upstream  boundary  condition  and  computations  were 
made  for  the  isolated  base  region  only.  This  allowed  a  large  number  of  grid 
points  to  be  used  in  the  base  region.  The  results  of  jet-on  calculations  for 
both  grids  are  shown  In  Figure  15.  Both  of  these  jet-on  calculations  give 
practically  the  same  base  pressure  variation  along  the  base.  At  the  nozzle 
exit  plane,  the  base  pressure  shown  is  the  specified  jet  pressure.  For  the 
remaining  portion  of  the  base,  the  effect  of  the  jet  has  been  to  reduce  the 
base  pressure.  The  small  oscillations  in  the  base  pressure  distribution  along 
the  base  are  located  near  the  two  corners:  one  at  the  base,  and  the  other 
formed  by  the  base  boundary  and  the  nozzle.  These  oscillations  are  not  physi¬ 
cal  and  are  believed  to  be  due  to  the  numerical  difficulties  associated  with 
the  boundary  conditions  at  these  corners.  The  base  pressure  distribution  in 
general  however,  is  a  realistic  representation  associated  with  the  base  region 
flow  field. 

Additional  qualitative  features  of  the  base  region  flow  field  are  shown  in 
Figures  16,  17  and  18  for  the  refined  grid.  Figure  16  shows  the  pressure  dis¬ 
tribution  along  the  surface  of  the  model  and  at  Z/D  =  .5  downstream  of  the 
base  corner  as  a  function  of  axial  position.  The  pressure  distribution  shows 
the  rapid  expansion  which  occurs  at  the  base  corner  and  the  recompression  that 
occurs  downstream  of  the  base.  Figure  17  is  a  pressure  contour  plot  in  the 
wake  region  and  clearly  shows  the  expansion  at  the  base  and  recompression 
shock  downstream  of  the  base.  Figure  18  shows  the  Mach  contours  in  the  near 
wake  flow  field.  The  strong  shear  layer  as  well  as  the  jet  boundary  in  the 
base  region  are  apparent.  It  also  seems  to  indicate  the  presence  of  a  Mach 
disk  at  X/D  »  9.75.  The  grid  spacings  in  that  vicinity  are  however,  not  as 
fine.  Thus,  a  sharp  picture  of  the  Mach  disk  is  not  seen. 


VII.  SUMMARY 

A  recently  developed  numerical  capability  has  been  used  to  compute  the 
missile  base  flow  in  the  presence  of  a  centered  propulsive  jet.  The  thin- 
layer,  compressible,  Navier-Stokes  equations  are  solved  using  a  time-depend¬ 
ent,  implicit  numerical  algorithm. 

Numerical  computations  have  been  made  to  predict  the  flow  field  over  a 
missile  afterbody  for  both  jet-off  and  jet-on  conditions.  The  free  stream 
Mach  number  is  1.343  and  ^he  jet  exit  Mach  number  is  2.7.  The  exhaust  jet 
static  pressure  is  2.15  times  the  free  stream  pressure  and  the  conical  nozzle 
half  angle  is  10°.  Qualitative  features  of  the  near  wake  flow  field  have  been 
shown  in  the  form  of  velocity  vector  and  contour  plots  of  stream  function, 
Mach  number  and  pressure.  Quantitatively,  the  effect  of  the  jet  has  been 
shown  to  reduce  the  base  pressure. 

The  experimental  results  are  not  available  pending  a  sealed  envelope  com¬ 
parison  with  computations  by  various  groups.  Comparison  with  the  experimental 
measurements  for  the  same  model  and  same  conditions  will  be  made  by  MICOM. 


Figure  1.  Schematic  Illustration  of  Flow  Field  Segmentation 
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LIST  OF  SYMBOLS 

a 

speed  of  sound 

a- 

free  stream  speed  of  sound 

• 

A 

cross  sectional  area  at  the  base 

CP 

specific  heat  at  constant  pressure 

J" 

- 

CP 

pressure  coefficient,  2(p  -  pj/p.u^ 

• 

D 

body  diameter  (57.15mm) 

e 

total  energy  per  unit  volume/p„a£ 

A 

q 

vector  of  dependent  variables 

•  . 

A  A 

E,  F 

flux  vector  of  transformed  Navier- Stokes  equations 

A 

H 

n-invariant  source  vector 

; 

I 

identity  matrix 

9 

J 

Jacobian  of  transformation 

M 

Mach  number 

M. 

free  stream  Mach  number 

•  . 

— 

P 

pressure/p^a* 

-  • 

P. 

free  stream  pressure 

Pr 

Prandtl  number,  p„c„/k 

9. 

R 

body  radius 

$•: 

Re 

Reynolds  number,  P^ajl/u.. 

S 

viscous  flux  vector 

• 

t 

physical  time 

•  .*.  % 

u,v,w 

Cartesian  velocity  components/a,,, 

V* 

u. 

free  stream  velocity 

• 

1 

u.v.w 

Contravariant  velocity  components/a,,. 

x,y,z 

physical  Cartesian  coordinates 

a 

angle  of  attack 

•  _ 
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LIST  OF  SYMBOLS  (continued) 


y  ratio  of  specific  heats 

k  coefficient  of  thermal  conductivity/^ 

coefficient  of  thermal  conductivity  at  free  stream  conditions 

u  coefficient  of  viscosity/u^ 

coefficient  of  viscosity  at  free  stream  conditions 

5,n,c  transformed  coordinates  in  axial,  circumferential  and  radial 
di recti  on s 

p  density/p,,, 

free  stream  density 

t  transformed  time 

<P  circumferential  angle 

Superscript 

*  critical  value 
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